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TECHNICAL NOTE 2677 


WING-BODY INTERFERENCE AT SUPERSONIC SPEEDS 
WITH AN APPLICATION TO COMBINATIONS 
WITH RECTANGULAR WINGS 
By Jack N. Nielsen and William C. Pitts 


SUMMARY 


An exact theoretical method is developed that permits the determina- 
tion of the pressure field of a wing-body combination having a circular 
body and a wing with supersonic leading and trailing edges. Detailed 
calculations have been performed for wing-body combinations composed of 
rectangular wings mounted at incidence on bodies at zero angle of attack 
for effective chord-radius ratios of 4 or less. For large effective 
chord-radius ratios some asymptotic results have been obtained. It was 
determined that for the family of combinations having an effective chord- 
radius ratio of 4 the area of the wing blanketed by the body does not 
generate any lift itself but rather acts to support the lift generated 
by the exposed wing, and that the body is less than 50 percent effective 
in reflecting lift back to the wing. For chord-radius ratios less than 
the relative amount of reflection increases. 

The significant fact was determined that for rectangular wing-body 
combinations for which the effective chord-radius ratio is greater than 4 , 
most of the loss of lift due to interference can be estimated from the 
first term of the Fourier series used in the analysis. This fact was 
used to determine asymptotic lift results for the region where no exact 
calculations were made. The asymptotic expressions, together with the 
calculations, allowed the construction of design charts showing the lift 
and center-of -pressure location of the exposed wing panels as a function 
of effective aspect ratio and effective chord-radius ratio. The charts 
show that as a result of interference the lift on the exposed wing panels 
in combination with the body can be reduced as much as 15 percent below 
the value for the wing panels joined together, and that the center of 
pressure of the exposed wing panels can move forward as much as percent 
of the wing chord. 


INTRODUCTION 


In recent years the problems of supersonic wing-body interference 
have occupied the attention of many workers in aerodynamics . The large 
amount of effort expended on the subject is a result of the important 
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effects that interference can have on the over-all aerodynamic char acter - 
istics of wing-body combinations. The trend, toward using large bodies 
and smal 1 wings at supersonic speeds , especially for missiles, is the 
prime reason for the increased importance of wing-body interference at 
these speeds. 

Much significant work has already been done in the field. In ref- 
erence 1, Spreiter has shown that, when a wing-body combination is slender 
in the sense of his paper, simple expressions for the lift and moment 
coefficients can be derived. These results were obtained by reducing a 
three-dimensional problem for the wave equation to a two-dimensional 
problem for Laplace's equation. Another approach is that of simplifying 
the differential equation by using conical boundaries . Following this 
approach, Browne, Friedman, and Hodes in reference 2 obtained a solution 
for the pressure field of a wing-body combination composed of a flat 
triangular wing and a cone both with a common apex. The use of all- 
conical boundaries reduces the problem to one of conical flow for which 
powerful methods of solution are available. 

Several investigators have presented methods for determining the 
pressure field, including the effect of interference, acting on wing-body 
combinations employing circular fuselages and wings not necessarily 
slender. In reference 3, Ferrari has given an approximate method of 
obtaining the "interference of the wing on the streamlined body, assuming 
that the induced field generated by the wing is that which would exist 
around the wing if it' were placed in the uniform stream alone . " Simi- 
larly, the Interference of the body on the wing has been determined. The 
results of Ferrari thus represent a first approximation and, while a 
second approximation using the method is possible in principle. It appears 
that too much labor would be involved. 

Another method for estimating the effect of interference on the aero- 
dynamic properties of wing-body combinations which are not necessarily 
slender is given in reference 4. In this reference the method is applied 
to dete rminin g the drag of symmetrical wing-body combinations; it is also 
applicable to the calculations of the lifting pressures acting on combi- 
nations employing wings with supersonic edges . In reference an essen- 
tially new method of solving a wide class of wing-body interference pro- 
blems has been presented. The method is based on decomposing the inter- 
ference of a wing-body combination into a number of Fourier components 
and solving the problem for each component in a manner similar to that 
used by von Ka r man and Moore in reference 6 for bodies of revolution. 

To su mmar ize the present situation. It can be said that the existing 
interference methods fall into two categories . The first category 
includes those attempts to solve difficult boundary-value problems as in 
references 3* 4, 5, and. 7- These methods have the shortco ming s that they 
all involve a great deal of labor, and that, with the exception of ref- 
erence 5, they represent approximations to the true linear-theory solu- 
tions. The second category includes simpler methods such as those of 
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references 1, 2, and 8 and existing approximate engineering methods. 

These methods have the disadvantages of being applicable only to special- 
ized configurations or of being based on assumptions of unknown validity. 

At the present time, there is a definite need for a simplified gen- 
eral theory of wing-body interference that will give results of engineer- 
ing accuracy for a wide range of wing aspect ratio, taper ratio, leading- 
edge sweepback angle, and span -diameter ratio. The reasons that none has 
yet been developed are twofold: first, the lack of systematic experi- 

mental results makes an empirical theory difficult; and second, incomplete 
knowledge of the mechanism of wing -body interference, because of the 
shortage of exact solutions, makes it difficult to develop a reliable 
rational method. The development of a simple general, theory of wing-body 
interference seems thus to be dependent on obtaining a few exact solutions 
that will give sufficient insight into the mechanism of interference to 
permit valid simplifying assumptions . Such solutions are also necessary 
to assess the validity of the assumptions underlying present theories. 

In reference 5> a solution has been obtained to a wing-body problem 
that is exact within the limitations of linear theory. It is the pur- 
pose of this paper to present a resume of reference 5 "with special 
emphasis on the physical principles underlying the interference. It is 
hoped that the insight gained into the mechanism of interference in this 
paper will be sufficient to allow the development of a simple, reliable 
engineering method. 

In this paper, a physical description of the interference problem 
will be given, and an Illustrative example of a rectangular wing at inci- 
dence on a round body at zero angle of attack will be presented' in detail. 
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SYMBOLS 

body radius 

aspect ratio of wing formed by joining exposed half -wings 
together 

chord of rectangular wing 
effective chord-radius ratio 

•minimum wave drag coefficient of wing based on exposed 
wing area 

chord at wing -body juhcture 
chord at. wing tip 
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f 2n (x) velocity amplitude function of n’th Fourier component 

g 2n (x) axial strength function of n’th Fourier component 

hg^x ) axial strength function for unit step in f 2Q (x) at 

origin 


K 


hie 

hi 


■^■gn( T j r ) 

hi 

hie 

L 

L nn 

m 

M 

M an( x ) ) w 2n( x ) 
n 


function used in determining P^ 

lift on exposed half -wings joined together 

lift on exposed half -wings in combination with body 

lift 

lift due to n'th Fourier component 
odd integer 

free -stream Mach number 
characteristic functions 
number of Fourier component 


P 

P pTl 

r^x 


s 

t 



pressure coefficient 

interference pressure coefficient due to n'th Fourier 
component 

free-stream dynamic pressure 

cylindrical coordinates: y = r cos 0, z = r sin 6 
(See fig. 1. ) 

semispan of wing-body combination 
maximum wing thickness 

thickness ratio of double-wedge wing 



axial, lateral, and vertical perturbation velocities, 
respectively 


u,v,w 
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c P2n*( I V 0 > x ) 

t(x,r) 


Cartesian coordinates: x, axial coor din ate; y, lateral 
coordinate; z, Vertical coordinate 
(See fig. 1.) 

center of pressure of wing lift measured from wing lead- 
ing edge 

free-stream velocity 
body angle of attack 
local angle of attack 

upwash angle of body-alone flow- 
wing angle of attack 

a/m 2 -1 

effective aspect ratio 
sin _1 x 

taper (s;) 

dummy variables of integration 
sweep angle of wing leading edge 
interference perturbation velocity potential 
combination perturbation velocity potential 
wing-alone perturbation velocity potential 
n’th Fourier component perturbation velocity potential 
fundamental solution for n'tb Fourier component 
undetermined function of x and r 

PHYSICAL PRINCIPLES 



Prior to a mathematical formulation of the wing-body interference 
problem, it is well to define interference and to explain bow it arises. 
With a stationary wing or a stationary body in a uniform parallel flow. 
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there are associated the wing-alone and body-alone flow fields. If the 
wing is imm ersed in the body-alone flow field, the fluid velocity due to 
the body-alone field will not, in general, be tangential to the wing (or 
conversely). For this reason, the sum of the body-alone plus wing-alone 
flow fields will not be the flow field for the body and wing together. 

The difference between the flow field of the body and wing together and 
the sum of the body-alone and wing-alone flow fields is defined to be the 
interference flow field. 


Effect of Forebody on Wing 


The effects of wing-body interference on the flow field of a wing- 
body combination are illustrated by considering separately the effects of 
each component on the others. For the purposes of this discussion 
figure 1 shows a wing-body combination divided into the part in front of 
the leading edge of the wing-body juncture, henceforth called the nose, 
the winged part and the part behind the wing trailing edge, henceforth 
called the afterbody. The wing may be twisted and cambered, but it is 
assumed to have supersonic leading edges. Consider now the flow as it 
progresses past the body- At the body nose the flow is that around a 
body of revolution, and it can be treated by existing methods such as 
those of references 1, 6, and 9* When the body is at angle of attack c&g, 
there is an upwash field in the horizontal plane of symmetry of the body. 
If the body is sufficiently slender, the flow field in a plane at right 
angles to the body axis corresponds to that around a circular cylinder 
in a uniform stream of velocity, V sin otg. This gives an upwash field 
in the horizontal plane of symmetry of the body of 

Ou = a B (l + a^y 2 ) (l) 

The effect of this upwash on the wing can be obtained by considering the 
wing to be at single of attack and twisted according to equation (l) and 
by applying the formulas of supersonic wing theory. The wing pressure 
field so obtained is exact, within the limitations of the theory, for 
that section of the wing outboard of the Mach line emanating from the 
leading edge of the wing-body juncture. If the wing is located close to 
the body nose so that there is a chordwise variation in the upwash field 
due to the body, then the wing is effectively cambered, and the solution 
is more difficult. However, for most wing-body combinations it is pos- 
sible to disregard the effect of the nose, and to assume that the wing • 
is attached to a circular cylinder that extends upstream indefinitely. 


Mutual Effects Between Body and Wing 


The mutual interference between the body and wing on the winged part 
of a combination causes an interference field acting on the body and on 
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the wing inboard of the Mach line emanating from the leading edge of the 
wing -body juncture. 

The wing-alone flow field does not, in general, produce flow tan- 
gential to the position to be occupied by the body surface. An inter- 
ference flow field must arise that cancels the velocity induced by the 
wing-alone flow field normal to the body while not changing the body 
shape. Alternately, the origin of the interference field can be explained 
in the following manner. The wing and body can be thought of as sources 
of pressure disturbances that radiate in all directions in downstream 
Mach cones. The wing disturbances which radiate toward the body are in 
part reflected back by the body onto the wing and in part transmitted 
onto the body giving rise to interference pressures. Likewise, the dis- 
turbances originating on the body pass onto the wing and affect-the pres- 
sures there. It is apparent that the determination of the interference 
pressure field on the body and on the wing inboard of the Mach line of 
the juncture is the crux of the wing-body interference problem. 


Mutual Effects Between Wing Panels 


To determine the region of influence of one wing panel on another, 
it is necessary to trace the path of a pulse ftom one wing panel across 
the body onto the other. The path traced across the body by the pulse 
originating at the leading edge of the wing-body juncture is the forward 
boundary of the region of influence of one wing panel on the body. (See 
fig. 1.) It is clearly the helix intersecting all parallel elements of 
the cylinder at the Mach angle. The boundary crosses the top of the body 

a distance of — */ M 2 -l downstream and reaches the opposite wing-body 

juncture a distance Jta a/M 2 -1 downstream. A pulse originating at a 
point on one wing panel and traveling to a point on the other panel can 
travel around the body on its surface to the opposite juncture and then 
along the wing to a given point, or.it can leave the body tangentially 
before reaching the opposite wing juncture in a straight path to the 
point. The second means of transmitting the impulse is shorter in dis- 
tance than the first and is the one which deter min es the forward bound- 
ary of the region of influence of one wing panel on the other. Applying 
this consideration to the pulse originating at the leading edge of one 
wing-body juncture, it is easy to show that the f orward boundary of the 
region of inf luence of one wing panel on the opposite wing panel is 
given by the equation 




cos" 1 



a + 




( 2 ) 


This bo undar y is also shown in figure 1, and it becomes parallel to the 
Mach direction at distances far from the body. 
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Effects on the Afterbody 


As far as the interference effect of the body on the wing is con- 
cerned, it is confined to the winged part of the combination, but the 
effect of the wing on the body is felt also on the afterbody. If the 
downwash were known everywhere in the wing wake, then the wake could be 
considered as an extension of the wing with twist and camber. The wing 
wake and afterbody could then be incorporated with the winged part of 
the combination and treated in the same manner. However, the actual 
downwash pattern in the wing wake depends on the interference effect of 
the body on the wing. It is thus apparent that the solution of the after- 
body problem requires that the interference problem for the winged part 
of the combination be solved first. Only the winged part of the combina- 
tion is analyzed in detail in this report. 

MATHEMATICAL FORMULATION OF THE PROBLEM 


Throughout the analysis, the body radius is taken as unity and M 2 
is taken as 2 so that p = 1. Any formula can be generalized to any body 
radius by dividing all length symbols by a, and to any Mach number by 
dividing all streamwise lengths by p, by multiplying al 1 pressure and 
lift coefficients by p, and leaving all potentials, lift forces, and 
span loading unaltered. It is necessary to specify the wing alone before 
any detailed interference calculation can be carried out. However, in 
the theoretical solution of the problem the wing-alone definition 1 b 
arbitrary. The flow field about the combination does not depend on the 
definition of the wing alone. 


The analysis is confined to the cases for which the aspect ratios 
are sufficiently large that the Mach lines from the leading edges of the 
wing tips do not intersect the wing -body juncture. This is the case if 
the aspect ratio obeys the following inequality: 


PA^: 


4 


1 + X 


)C 


tanA 


+ 1 


) 


Under such circumstances, the tips have no influence on the wing-body 
interference of the winged part of the combination. 


The wing-body combinations being considered in this report are those 
possessing a vertical plane of symmetry. To eliminate any questions of 
the symmetry of the boundary conditions above and below the wing, the 
assumption is made that the wings have supersonic leading and trailing 
edges. This assumption, together with that concerning tip effects, is 
sufficient to insure the fact that pressure communication between the 
upper and lower wing surfaces can have no effect on the interference on 
the winged part of the combination. Thus it is immaterial whether the 
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problem at hand is symmetrical or antisymmetrical with regard to the 
horizontal plane of symmetry. Attention will henceforth be focused above 
the horizontal plane of symmetry. 

The assumptions with regard to tip effects and the sweep of leading 
and trailing edges enable the drag of a symmetrical wing-body combination 
at zero angle of attack to be obtained from solutions for lifting sur- 
faces. Suppose the problem solved is that for a flat wing at angle of 
attack oty mounted on a body at zero angle of attack. The pressures 
are antisymmetric with respect to the horizontal plane of symmetry. How- 
ever, considering the same pressures as positive and symmetrical yields 
the solution for the pressure field of a single-wedge wing of half 
angle oty at zero angle of attack. This fact is of use in obtaining the 
drag coefficient at zero angle of attack, as will be illustrated. 1 


General Decomposition of Boundary-Value Problem 


The general case of the wing and body at different angles of attack 
is considered. Following the suggestion of Lagerstrom and Van Dyke in 
reference 10, the problem can be broken down into several simpler prob- 
lems on the assumption that the wing boundary conditions can be applied 
in the plane of the wing and those for the body can be applied on a 
cylindrical surface. Figure 2 shows the decomposition of the wing-body 
combination into three components. Component (a) is simply the body 
alone, which creates an upwash field in the region that is to be 

occupied by the wing. Components (b) and (c) are combinations with wings 
of the same plan form; but while component (b) has the wing-alone angle 
of attack, component (c) has a wing with angle of attack -a^. The sig- 
nificance of this particular method of decomposing the general wing-body 
problem is that component (a), the body alone, can be solved by known 
methods and components (b) and (c) with bodies at zero angle of attack 
can be solved by the methods of this report. Henceforth the analysis is 
confined to combinations with bodies at zero angle of attack. 

Wing alone and interference potentials .- Consider now a combi n ation 
with the body at zero angle of attack and let qu be its potential. 

(See fig. 3 . ) This potential can be considered the sum of a wi ng -alone 
potential qfy and of an interference potential qp. 

9 C = \ + 9 ^ 

The essential problem is to determine cp. First, select a convenient 
way of extending the wing through the body to form the wing alone, thereby 
specifying 9^. The wing-alone flow field in general produces veloc- 

ities — normal to the surface that will enclose the circular cylinder 
or 

as illustrated in figure 3 for the region above the wing. In figure 3 
and the subsequent figure, all bodies are shown as cylinders parallel 
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to the x axis. While the bodies of the actual configurations in same 
cases are slightly distorted cylinders, they are nevertheless shown as 
true cylinders . This procedure is compatible with the fact that the 
boundary conditions are to he applied on a true cylinder. The value 

of varies with 0 and with x. This means that a body conforming 

dr 

to the wing-alone flow field is distorted in a complicated fashion. Now 
since the body must be circular, there must arise an interference poten- 
tial cp that identically cancels at the body surface, thereby 

or 

straightening it. 


^=-4 at r=l (4) 

or or 


There are two other conditions to be fu lf illed by T* It must not dis- 
tort the shape of the wing when added to cpy to produce Tq* Thus 
when 0 = 0 , 


dT _ 

d0 


(5) 


or cty = 0 for the interference combinations as Bhown in figure 3 * The 
last condition is that the interference potential must be zero ahead of 
the winged part of the combination. 

cp = 0 x5-0 (6) 

Equations (4), ( 5 ), and ( 6 ) are the essential boundary conditions on T. 

The normal velocity to be induced at the body surface by the 

dr 

interference potential can be analyzed at any given streamwise position 
as a Fourier cosine series. The amplitudes of the various Fourier 
cosine terms, f 2 Q (x), vary with x, the streamwise distance. Thus, 


^ f 2 n( x ) cos 2110 = - ^ at r = 1 (7) 

n=o 

Only even multiples of 0 are considered because of the vertical plane 
of symmetry. Consider that the interference potential is decomposed into 
a series of potentials such that each cancels one Fourier component of 
the velocity at the body surface, that is, 

00 

9 - Z ^ 

n=o 


( 8 ) 
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with 


dq>an 

Sr 


= f 2 n(x) cos 2n 9 


at r = 1 


(9) 


Then the combination giving the interference potential $ can be decom- 
posed in a series of combinations each giving one of the 9 pn values. 
The decomposition is illustrated in figure For n = 0, 


5( P 0 

Sr 


fo( x ) 


and there is no variation of the normal velocity, pressure, or potential 
with 0. Thus the first interference combination is a body of revolution. 
The pressure field acting on the body of such a combination can be deter- 
mined by the method of reference 11. This n = 0 interference combina- 
tion has the very simple significance that its flow normal to the 

r = 1 cylinder, subtracted from EEE reduces the flow across the 

Sr Sr 


body to zero when averaged from 0 = 0 to 9 = n at any streamwise loca- 
tion. For n = 1, 


Sqpg 

Sr 


f 2 (x) 


COB 20 


and the normal velocity, pressure, and potential will vary as cos 20. 

To summarize briefly, it has been shown that the general inter- 
ference problem of a body and wing at different angles of attack can be 
broken down into wing-body problems with bodies at zero angle of attack 
as shown in figure 2. Combinations with the body at zero angle of attack 
are decomposed into wings alone plus interference combinations as in 
figure 3- The interference combinations are finally decomposed into 
their Fourier components as in figure 4. 

A general method for determining the characteristics of any Fourier 
component will now be given. It will be shown that good accuracy can be 
obtained for the interference potential with few Fourier components. 


Method of Solution 


The problem to be solved is that of a supersonic wing-body combina- 
tion subject to the conditions already mentioned, but with the wing and 
body possibly at different angles of attack. This problem is reduced to 
a body-alone problem and two wing-body probl e ms with the body at zero 
angle of attack as shown in figure 2., The body-alone problem can be 
solved by existing methods such as references 6 and 9* The procedure 
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necessary to solve either wing-body problem as given in reference 5 is 
now summarized. 

Wing-alone potential and velocity amplitudes .- The first step in 
the solution is to deter min e the wing-alone potential, cp^j- The methods 
of wing theory available for doing this will not be discussed here. 

From the value 0 of the value of the normal velocity induced at 

the r = 1 cylinder follows readily. This normal velocity distribution 
is expanded in a Fourier cosine series of even multiples of 0. 

CO 

^ f 2n (x) cos 2n0 at r = 1 (10) 

n=o 

The functions fpn(x) are called the velocity amplitude functions and 
are part of the boundary conditions for the interference potential. 

Interference potential .- The solution for the interference potential 
is obtained by summing the potentials for the series of Fourier component 
potentials. The determination of the potential for any one Fourier com- 
ponent (as illustrated by fig. k) has been carried out formally in ref- 
erence 5 by Laplace transform theory. It was shown that the formal 
mathematical solution is equivalent to distributing fundamental solutions 
for cp along the axis of the combination in a manner similar to that of 
von Karman and Moore in reference 6. Let cp 2n *(x-£,r,0) be the funda- 
mental solution located at x = | to be used in determining the poten- 
tial of the n'th Fourier interference component. Let g 2n *(i) be the 
distribution of this solution along the body axis from x = -1 to x = oo 
to satisfy the boundary conditions. Then the potential of the n'th 
Fourier component is 



c P 2 n( x >r,0) = jT 9an*( x - fi > r >0) S2n*^) 

First, the set of fundamental solutions cp 2n *(x,r,0) will be presented, 
and then the method of obtaining g 2 n*(S) be discussed. 

Axial strength function .- The boundary conditions to be satisfied 
by each Fourier component of the interference potential are given by 
equations (5), (6), and ( 9 ). The first two boundary conditions are 
satisfied by choosing an appropriate set of fundamental solutions. Such 
a set of fundamental solutions, derived in appendix A, is 


. ^ 2 n *( x “^r,0) 


cos 2n0 d 211 " 1 

r 2H rbr 2 * 1-3 - 


[(x-|) 2 -r a ] 



( 11 ) 


The fundamental solution is located at the point x = £, and the value 
of l varies from -1 to +». For n = 0 an integration is to be per- 
formed. The last boundary condition is satisfied by distributing 
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fundamental solutions, as given in equation (ll), along tie body axis in 
variable strength. In appendix B tbe following equation is derived for 
dete rmining tbe axial strength function 

SaiW = M *n(*-5) « (12) 

The evaluation of g^fe) from this equation is done numerically or 
graphically. The values of f^^x) are obtained on the basis of equa- 
tion (10). A table of the functions M 2n (x) are included as table I for 
use in equation (l 2), and figure 5 illustrates these functions. The 
properties of M^^x) are discussed in reference 5« 

* 

Potential or span loading .- Once the axial strength distribution of 
the fundamental solution is known, the determination of the interference 
potential or pressure follows readily. Consider the general point P, 
shown in figure 6, where the potential is to be determined. The funda- 
mental solution along the axis from -1 to x-r influ ences the potential 
at P which is given as 

tjx.r.e) = f X ’ T [(x-6) 2 -r=:r-i g2n «( s ) dj 

231 J_ i r® dx 31-1 . 

(13) 

By making a unit translation of the variable of integration so that the 
lower limit is zero, there is easily obtained 


<P 2n( X , r , 9 ) 


cos 2n9 



x-r+i 


a 211 - 1 

dx 2n_1 


[(x-i+i)2-r 2 f“ s &*({) at ' UM 

Equation (l4) serves to determine the potential anywhere in the flow 
field once the axial strength function g 2n (x) has been obtained. Since 
the integrated lift per unit span up to a given streamwise position x 
is proportional to cp at that position, equation (l4) also gives the 
span loading. 


Pressure coefficient .- On the basis of linear theory, the pressure 
coefficient is 


P = 


2u _ _ 2 d<P/Sx 
V 


V 


(15) 
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With the aid of equation (15) it is readily shown from equation (l^) that 


where- 


r sn 


= -2 


cos 2n0 


l 


x-r+i 


K 2 el ( t J r ) g^Cx-r+l-T) dT 


Kpri(T,r ) 


1 S £n [T(T + 2r)] 2n ~i 

r 211 ^ T 2 n 


(16) 


The procedure in obtaining the pressure coefficient is first to 
obtain the axial strength function from equation (12) by numerical or 
graphical integration, and then to obtain P pri for any point in the 
flow field from equation (l6) by the same means. 

Body pressure coefficient .- The determination of the pressure coef- 
ficient by the method outlined requires two numerical or graphical inte- 
grations. In appendix C a simpler method involving only one integration 
is derived for the body pressure coefficient. The result for the pres- 
sure coefficient is 


PgrL = y cos 2x19 f an( x ) “ | cos 2n0 J I )W 2n (x-& ) dl (17) 

at r = 1. The functions W ai (x) are universal functions. They are 
tabulated in table II, and are plotted in figure 7 • 

Combination potential .- The combination potential is obtained by 
simply ad din g the wing-alone potential to the interference potential cp. 
The interference potential cp is obtained as already described by add- 
ing together the interference potentials of a number of Fourier compo- 
nent combinations. The lift or pressure coefficients for the combina- 
tion are similarly obtained simply by adding to those for the wing alone 
those for the interference combination. 

Referring again to figure 2, the detailed method for solving either 
cases (b) or (c) has been summarized. To solve the general case of both 
body and wing at angle of attach, both (b) and (c) must be solved. The 
wing alone for each case is different. The simplest interference problem 
is that for the body at zero angle of attack. Under these circumstances , 
(% is zero, and only case (b) must be solved. This case will be calcu- 
lated for a rectangular wing-body combination as, an illustrative example. 
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PRESSURE FIELD ACTING ON A COMBINATION WITH RECTANGULAR WINGS 


One of the simplest wing-body combinations of technological impor- 
tance is that formed by a rectangular wing mounted at incidence on a cir- 
cular body at zero angle of attack. The complete pressure field acting 
on such a combination will now be determined. The results are inter- 
preted in terms of lift and center of pressure presented as a function 
of effective aspect ratio and effective chord-radius ratio. As pre- 
viously mentioned, the wing alone can be specified in any convenient 
manner and, for the purpose of the example, the wing alone is taken as 
the rectangular wing extending straight through the body from side to 
side. Although the analysis as carried out is for M = the results 

are presented in a form applicable to a range of Mach numbers. The steps 
in performing the calculation are: (l) to determine the wing-alone 

potential; (2) to determine the velocity amplitude functions, f^Cx); 

(3) to determine the axial strength functions, g^(x); and (4) to deter- 
mine the potential or pressure, as desired, anywhere in the field. No 
tip effects are considered until the results are presented as a function 
of wing aspect ratio. 


Wing-Alone Potential 


The wing-alone flow, exclusive of tip effects, can be determined 
from the Ackeret theory. The flow at a spanwise station out of the 
region of influence of the wing tips is illustrated in figure 8. The 
potential for the flow above the wing is 


= Vx when z>x 


(18) 


— H. is uniformly -a^V. 


= Vx + a^vCx-z) when z<x 

^Pw 

dy 

The downwash causes a 


(19) 


is zero and the downwash 


flow normal to the sur- 


The sidewash produced by such a potential 
8z 

face r = 1 in amount -cfyV sin 9 . This means that for a body conform- 
ing to the wing-alone flow, the deformation is zero at the wing-body 
junctures and a maximum on the top of the body. The interference combi- 
nations when added to the deformed body straighten it out Fourier com- 
ponent by Fourier component. 


Fourier Amplitudes of Body Normal Velocity 


The Fourier amplitudes of the normal velocity induced by the wing- 

alone potential at the body is determined by expanding — - at r = 1 

Sr 
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in a Fourier cosine series of even multiples of Q. The normal velocity 
distribution is shown in figure 9* For x>l the body is totally 
immersed in the wing downwash field. With the usual equation for obtain- 
ing the Fourier amplitudes of a function, there is obtained 


/->sin -1 x 

f Q (x) = - / OgV sin 6 d9 

'-'n 


fpn ( x ) 


The integrations give 


If 

*1 


sin _1 x 


OgV sin 6 cos 2n0 d0 


f Q (x) = ^1 - l-x 2 ^ when x<l 


f Q (x) = when x>l 


( 20 ) 

( 21 ) 


( 22 ) 

(23) 


- St 


2 cos (2n-l)S _ 2 cos (2n+l)5 _ 4 

2n-l 2n+l kn 2 -l 


when x < 1 


Fpn( x ) 


-Wag 

rt(4n 2 -l) 


when x > 1 


(24) 

(25) 


where 5 = sin -1 x. The ,^2n( x ) functions are shown in figure 10. The 
constant values of f 2n( x ) for x>l are noteworthy. 


Axial Strength Functions 


The calculation of the axial strength distributions has been carried 
out with the aid of equation (12). The integrations were performed 
numerically by means of integration formulas given in reference 12. Two 
separate cases occur, x<l and x>l. For x<l a typical plot of the 
integrand and its components is shown in figure 11(a). "It will be 
observed that the Mq(x-£ ) function has a singularity at the upper limit 
of integration. Since the singularity is of the square root type, it is 
easily handled. For x>l the determination of go(x) is complicated 
by the fact that f 0 ’(i) has a square root singularity at 1=1, but 
simplified by the fact that f 0 '(i) is zero for i >1. A typical plot 
of the integrand and its components for x>l is shown in figure 11(b). 
The singularity of M Q (x-i) has been replaced by that of f 0 ’(|) at the 
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upper limit of integration. For x = 1 the square root singularities 
of Mq(x-§ ) and f 0 '(£) both occur at 1=1, reinforcing each other in 
such a manner that g^x) has a logarithmic singularity at x = 1. In 
fact, it can be shown that 

( ' 1)n (*) 

in the neighborhood of x = 1. The function g^Cx) minus the logarith- 
mic term passes smoothly through x = 1. 

A plot of the g 2 n(x) functions for n = 0 to n = 3 is shown in 
figure 12. The main effect noted is that the number of times the curve 
crosses the zero axis increases as the order of the g^n(x) functions 
increases. The singularities at x = 1 are also illustrated. 


Pressure Distributions 


Obtaining the strength functions of g 2n (x) is tantamount to getting 
the potential or pressure coefficient anywhere in the flow field as given 
by equation (l4) or (l6). If only the span load distribution is required, 
then one calculation using equation (l4) suffices for each spanwiBe 
station; but if the detailed pressure distribution is required, then 
many calculations using equation (l6) must be made for each station. The 
detailed pressure distributions have been calculated using equations (l6) 
and (17). In equation (l6), the factor 

WT,r) = -ij ^ [t(t^t)] 211 ^ 

appears inside the integral as an influence coefficient for converting 
axial strength function into pressure coefficient. The factor plays an 
important role in determining the accuracy of the numerical work. In 
figure 13 plots of Kq ( r, r ) and K 2 (-r,r) are presented showing the singu- 
larities of KgjjCTjr) at the origin of r. These singularities, together 
with those of g Pn (x-r+l-T) at t = x-r, bring about two different cases 
in the determination of the pressure coefficient. If x<r, only 
the K 2 n (T,r) singularity occurs in the integration as illustrated by 
figure l4(a) . However, for x>r there are two singularities as shown 
in figure l4(b). For x = r the two singularities come into confluence, 
but they still produce a finite pressure coefficient. 

Interference pressure distribution .- The interference pressure dis- 
tributions have been calculated for the first four Fourier components 
and are presented in figure 15 . In this figure the abscissa is propor- 
tional to distance behind the Mach line originating at the leading edge 
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of tlie juncture , as Ill ustrated in the figure. Although the calcula- 
tions have been carried out for M = •/ 2, that is, p = 1, and for unit 
radius, they are generalized to all Mach numbers and body radii by replac- 

X X 

ing x-r+1 by — — — + 1 and P^n by pP^ as has been done in the 

figure. From the figure it is apparent that the cusps in the pressure 
distributions are propagated downstream along lines of constant 

X X 

i-l or x-pr; that is, along the downstream characteristics. As 

pa 8. 

the pressure distributions move outward from the body along the down- 
stream characteristics, they are distorted and decreased in magnitude. 

Increasing the order of the Fourier harmonics causes two important 
effects: first, the number of points of zero pressure is increased and, 

second, the pressure coefficient damps more rapidly. As a result of the 
first effect the contributions of the higher harmonics to the combination 
span loading are proportionally less than their contributions to the 
pressure coefficient; while as a result of the second effect the more 
remote a point is from the leading edge of the wing-body juncture, the 
fewer the number of Fourier components that must be included to obtain 
its pressure coefficient accurately. All interference pressure distribu- 
tions exhibit discontinuities in slope at — — £. + 1 = 1. This behavior 

Ba a 

is a consequence of the fact that the body becomes totally immersed in 
the wing-alone flow field for this condition. When the pressure distri- 
butions of the various Fourier components are added together to obtain 
the interference pressure distributions, the discontinuities in slope 
tend to cancel so that the pressure distribution for the combination will 
be smooth. 

A detailed examination of the interference pressure distribution 
for the first Fourier component illustrates several points of interest. 

The importance of the component arises from the fact that it accounts 
for most of the effect of interference on the span loading. The reason 
for this is that the pressure coefficients for n = 0 are of invariable 
sign. The effect of the first Fourier component is to reduce the veloc- 
ity induced normal to the body by the wing-alone flow field to zero 
average around the body for 0 = 0 to 0 = it at any streamwise location. 

For purposes of comparison with the exact results for n = 0, some 
approximate results have been included in figure 15(a). For values 

of p^-< 1 on the body, the Ackeret value of P Q (twice the local stream 

angle divided by p) is a close approximation to the true pressure coef- 
ficient. This is the result of the facts that the part of the body 
affecting the interference is effectively plane for points near the lead- 
ing edge of the wing-body juncture and that there is no variation of any 
quantities with 9 so that an approximate two-dimensional situation pre- 

Y 

vails. As — increases beyond unity on the body, there is a rapid 

P a 

decrease in the pressure coefficient below the Ackeret value due to the 
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effect of all disturbances In front of the point in question as repre- 
sented by the integral of equation ( 17 )* 


In reference 5> the following approximate results were obtained for 


small and large values of 




\Pa 


/ 


for the pressure coefficient: 




2n 


cos 2nd 

a ,/r" 


#--£ + 1 

Pa a 


> t"£ + 1 ^ 0 


( 26 ) 



itx/pa'’ 


2L _ H + !_> a, 

pa a 


(27) 


Figure 15(a) shows that for values of ~ — ^ + 1<0.6 equation (26) is 
a good approximation for n = 0 although it is of little value for 
higher— order harmonics. There is a general tendency of P Q to approach a 
uniform value independent of r as + 1 becomes large, as shown 

by equation (27). The damping in the characteristic direction, although 
initially inversely proportional to the square root of r, is ultimately 
independent of r. 


The pressure distributions of figure 15 were determined by equa- 
tion ( 17 ) for r = 1 , and by equations ( 12 ) and ( 16 ) for r = 1 . 5 * 2 . 0 , 
and 3.0. However, the values for r = 1 were checked by equations (12) 
and ( 16 ) before these equations were applied to the higher values of r. 
For r = 1 equation ( 17 ) gives more accurate results than equations (12) 
and ( 16 ) since it requires only one numerical integration. The accuracy 

of both methods of calculation decreases aB n and — — + 1 increase 

8a a 

because the small values of the pressure coefficients under these condi- 
tions are the result of large counterbalancing influences. 


Pressure distribution in juncture of wing-body combination. - By 
adding the interference pressure coefficients of the various Fourier 
components to that for the wing alone, the pressure distribution for the 
combination is obtained. The addition has been carried out for the wing- 
body juncture using four Fourier components and, six Fourier components, 
and the results are presented in figure 16. The pressure coefficient 
with interference is less in magnitude than two, the value without inter- 
ference, showing that significant losses of lift occur in the wing-body 
juncture. A comparison of the results for four components and six com- 
ponents shows that four components give good over-all accuracy for all 


values of 


greater than 1. For small values of ■=— in the wing- 


x 

“ — Kicauci u iin.11 jl • 1!ujl PJiiau. vaiuco ujl 

pa Pa 

body juncture, the curvature of the body insofar as the flow is concerned 
is not large so that the body is effectively a vertical boundary on which 
a given distribution of normal velocity is producing an interference 
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field. Supersonic wing theory applied to this condition gives for the 
net interference pressure coefficient (reference 5) 


^ - when — > 0 ( 28 ) 

cfy 3rt^a pa 


It is clear that the calculated results can be joined smoothly to this 
result. Using the result of equation (28) enables satisfactory results 
to be obtained with four Fourier components . 


The critical region in the convergence of the solution is that near 
the leading edge of the wing-body juncture. The higher harmonics have 
their most important effect near here, and rapidly damp downstream along 
the body. Hence more and more Fourier components would be required to 

get accuracy for smaller and smaller values of — . However, with the 

pa 

result of equation (28), this extra work is unnecessary. 


One point of interest in figure 16 is the fact that when 

p& 

approximately 3 , the pressure coefficient increases in magnitude . 
is due to the fact that for 

pa 

wing is felt in the wing-body juncture 


equals 
This 

the influence of the opposite half- 


Pressure distribution on top meridian of wing-body combination . - 
The pressure distribution on the top meridian of the wing-body combina- 
tion is obtained in the same fashion as that at the wing-body juncture, 
the difference being -that the pressures due to the even number Fourier 
components have the same sign at the meridian as at the juncture, whereaB 
the odd numbered components have reversed signs. The pressure distribu- 
tions based on four and six Fourier components are shown in figure 17 . 


Several interesting effects are exhibited by the results. The step 

in the wing-alone pressure at — = 1 is effectively canceled by the 

~~ pa „ x / 

interference pressures of the Fourier components from = 1 to — = jt/2 

and for ~>n/2 the pressure increases rapidly and tends toward the two 

dimensional value. The effect of the interference pressure in canceling 

the effect of the wing alone on the top of the body from ~ = 1 


to 


pa 


— - = rt/2 is to be expected since the wing of the combination can 

pB 

have no effect on the top of the body unless — >tf/2, as has been 

pa 

already pointed out. If an infinite number of Fourier components had 
been taken, the pressure coefficients would be identically zero 
from — = 0 to ■—- = The general behavior in this regard is evi- 

pB pB 

dence of the plausibility of the calculated results. 
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The tendency of the pressures to approach an asymptotic value is also 
illustrated by figure 17- This asymptotic value represented by the sum 
of the wing-alone pressure plus the asymptotic results for the first 
Fourier component is given by the following equation: 

SJ -2 + —4— when — >2.4 ( 29 ) 

tfx/pa pa 

For ~>2.4. the results of this equation are in good agreement 
pa 

with the results of figure 17- 

Some evidence is furnished from the pressure calculations for the 
juncture and top of the body concerning the number of Fourier components 
necessary for accuracy. Comparisons made in figures 16 and 17 show that 
about four components are sufficient, and that the addition of two more 
is not worth the extra work. 

Pressure distributions for wing of wing-body combination .- The dis- 
tribution of the pressure acting on the wing of the combination can be 
determined in a manner similar to that for the wing-body juncture by 
adding to the wing-alone pressure those due to the Fourier components. 

The resultant pressure distribution for the wing based on four Fourier 
components is shown in figure 18. A drawing of the pressure field for 
the complete configuration is included as figure 19- For small values 
of — the higher-order oscillations in the pressure coefficient as 

shown in figure 1 6 have been ignored and the curves have been faired 
through them. 

Since the region of influence of the body on the wing is confined to 

the wing region downstream of the Mach lines emanating from the leading 

edge of the juncture, in front of this line the pressures are uniform 

at the two-dimensional value and behind the line there is a decrease in 

the magnitude of the pressure coefficient. If the body were a perfect 

reflector, that is, a vertical wall of infinite extent, then there would 

be no pressure loss. However, the pressure pulses originating on the 

wing are only in part reflected by the circular body. The efficiency of 

the body as a reflector is discussed subsequently in connection with span 

loading. The tendency of the pressure to increase in magnitude near the 

inboard trailing edge is due to the effect of the opposite wing panel 

which at the wing juncture is felt downstream of the point — = rt . 

pa 


Span Loading 


The span load distributions for a range of rectangular wing-body 
combinations with the body at zero angle of attack can be determined 
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from the pressure distributions of figure 18. For purposes of specify- 
ing the span loading the following equation is used: 



} (30) 



The span loading is taken as the quantity inside the brackets „ 


The pressure results of figure 18 are for values of the effective 

chord-radius ratio of 4 or less and for values of the effective aspect 

ratio of 2 or greater. Span loadings for any combination of (or c*) 

pa 

and pA in these radges can be obtained by integrating the pressure dis- 
tributions. The span loading evaluations have been made for e* = 4 
and pA>2. First the span loadings due to the various Fourier components 
are discussed, and then the span loadings for the actual wing-body combi- 
nations are presented. 

Fourier component cnmhl nations .- In figure 20, the contributions to 
the span loading for the first three Fourier components are shown. 

For n= 0 the pressure field does not depend on 0, being axially sym- 
metric, and a constant loading exists on the body. However, on the wing 
as the spanwise distance increases there is a decrease in the span load- 
ing due primarily to decrease in the length of chord over which the 
interference pressures act. The span loading due to the first Fourier 
component causes a loss of lift everywhere along the span. 


A comparison of the results of figure 20 for n = 0 and n = 1 shows 
that the first Fourier component accounts almost entirely for the effect 
of interference on the span loading of the combination. For the body 
this fact is even more true than for the wing. This fact is of consider- 
able importance since it gives a simple means of expen din g the lift and 
moment results to larger values of than those for which the pressure 

distributions have been calculated. Also, it suggests a simple means of 
minimizing the adverse effects of interference on lift as will subse- 
quently be pointed out. 


With the techniques of Laplace transform theory, it is possible to 
obtain asymptotic formulas for the span loadings of the various Fourier 
components. For the first Fourier component the following asymptotic 
result has been obtained by the standard methods of Laplace transform 
theory. (See reference 5* ) 


-2 



Po 

“W 


dx- 




+ — , — when — 
x/pa Pa 


— > 00 


(31) 
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The asymptotic result for the span loading given by this equation, when 
compared, with the results of the exact calculations in figure 20, is 
seen to be slightly high. However, for values of ^ greater than h 

the difference between the results decreases, and equation ( 31 ) thus pro- 
vides a satisfactory means of extrapolating the results of the present 
calculations for span loading to larger values of 


The asymptotic result has also been determined for the higher-order 
Fourier components as a matter of interest. The span loading is 



P pn , ^ 8 cos 2n0 

-= - dx-> 


, 2 n 


16 ( r 211 + ) (^n-l)Icos 2 n 0 


/ „ \4n 

m(2ni ) 2 (4n 2 -l)2 4 ' n “ 1 ( ~ ) 

\ P 3 -/ 


as 


pa ^ 


00 


(32) 


A comparison between the results of equation (32) and the exact solution 
for n = 1 in figure 20 shows agreement and corroborates the fact that 
the span loading of all but the first Fourier component is negligible 
for ^->4. It is also to be noted that the contribution to the loading 
of the first component given by equation ( 31 ) increases without limit 
as x — > 00 , whereas the span loading of the higher-order components is 
finite. 


Family of complete combinations .- To obtain the span loading for the 

family of combinations for which -S_ = 4, it is necessary to consider the 

pa 

loading of both the wing alone and the Fourier components. The necessary 
calculations have been carried out and the span loadings for the family 
of combinations based on one and four Fourier components are both shown 
in figure 21. The loading due to the wing alone is also shown. No 
effect of wing tips has been included. It is to be noted in figure 21 
that, whereas the loading on the wing due to its own pressure field is 
constant, there is some loss on the body because of the fact that the 
pressure field of the wing alone acts on the body only if x>pa sin 0. 
However, if an afterbody is included, some of the lift lost can be recov- 
ered. As has already been pointed out, the pressures due to the first 
Fourier components are positive on the upper half of the wing-body combi- 
nation and produce a loss of lift, as figure 21 shows. When the effects 
of four Fourier components are taken into account, the net lift is 
slightly higher than that for one Fourier component, but the difference 
is not significant. For most engineering purposes, one Fourier component 
is sufficient for determining the span loading when |£-> 4. 

Some insight into the mechanism of wing-body interference can be 
gained by comparing the span loading for the combination with those for 
two reference loadings: (l) the complete reflection case for which the 

blanketed area of the wing acts effectively at cty, and (2) the 
no-reflection case for which the blanketed area of the wing is supposed 



2b 


MCA TN 2677 


to act effectively at zero angle of attack. The span loading correspond- 
ing to the first case of complete reflection of the wing pressure pulses 
by the body is, in fact, the span loading marked "wing alone" in fig- 
ure 21. A comparison of this curve with that based on one or four 
Fourier components shows that the loading given on the assumption that 
the wing blanketed area is fully effective in lift is too optimistic. 
Under the conditions of the second reference loading, the sole purpose of 
the blanketed area is to support lift generated by the wings. A compari- 
son of the span loading for this case with the true loading shows that 
the average load on the body is well predicted, but that the loading on 
the wing is, of course, underestimated. A comparison of the true load- 
ing with those for the two reference cases reveals the interesting fact 
that the body is somewhat less than 50 percent effective in reflection 
for this particular family of configurations. 


Lift 


For values of — <4 the pressure distributions already presented 
pa 

are sufficient for obtaining span loading or lift on either the wing or 
body for al 1 combinations having sufficiently large aspect ratios to 
avoid effects of the tips on the wing-body interference. This is the 
case for pA>2. For values of larger than four, the lift results 

are presented in terms of a nondimensional parameter k w , defined as the 
ratio of the lift on the exposed half -wings in combination (exclusive of 
that on the body) to that on the exposed half -wings joined together. 


*w = 


% 


(33) 


For — >4 the value of k„- can be obtained by using the asymptotic 
Pa r \ 

form of the span loading given by equation ( 31 )* 


k w 





2 _ 1 
*c* 2c* 2 






— ~ — log 

itc^pA 



C* oo 


(34) 


The values of k v have been determined from the pressure distribu- 
tions of figure 18 for values of -£-<4 and from equation (34) for 

pa 

values of — >4. The effect of the wing tips has been taken into 
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account by utilizing reference 13. The results are shown in figure 22 

wherein kw is given as a function of for various effective aspect 

__ P®' 

ratios of 2 and greater. It should be borne in mind that the results of 
the figure are for a combination of body and rectangular wing or an all- 
movable, rectangular control surface with no gap. It is noted in the 
figure that the exact results for can be faired into the asymptotic 

results for thereby providing a design chart for engineerin g pur- 

poses for the entire range of p^. The curves of figure 22 illustrate the 
decrease of kw . as ^ increases at constant effective aspect ratio, 

and the slow increase of kw toward unity as the wing chord becomes 
very large. The loss of lift is most serious for f3A = 2, being about 
15 percent in the worst case. 

A practical point in connection with the loss of lift on the wing 
due to interference is that this loss occurs no matter what the body 
angle of attack, even though the calculations are made for Og = 0. It 
occurs either in the case of a wing mounted on a body or in the case of 
a deflected all -movable control surface. For wings with swept lead i ng 
edges for which all of the wing area lies in the region affected by the 
interference, even larger losses than occur with rectangular wings are to 
be anticipated. However, the loss of lift at the design condition can, 
at least in principle, be largely prevented by designing the fuselage 
so that it conforms to the first Fourier component in the wing-alone 
flow. This would involve contracting the fuselage above the horizontal 
plane of symmetry in a rotationally symmetric fashion and expanding a 
like amount beneath the horizontal plane of symmetry. Whether or not 
such a change would improve the lift-drag ratio can best be determined 
by experiment. 


Center of Pressure 


The center-of -pressure locations have been calculated for the same 
condition as the lift results of figure 22. The center-of -pres sure loca- 
tion in chord lengths beh in d the leading edge are presented in figure 23 * 
For large values of c*, an asymptotic result has been calculated for 
x c .p./c using the methods of Laplace transform theory and considering 
only one Fourier component . 


1 

x c.p. 2 

± — . <v _ 

C 


1 

3PA 


log 2c* 2 



as c* — ^ oo 


(35) 
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The values of X — ' *P .l. have been determined from the pressure distri- 
c 

butions of figure 18 for values of and by equation (35) for 

n \ P®' 

values of — >4. The loss of lift near the tips has been taken into 
Pa 

consideration. The exact results for — <4 have been faired into the 

pa 

asymptotic results for large values of ^ by dashed curves to provide 

an engineering design chart covering the entire range of ■£-. It is 

pa 

again mentioned that this chart is applicable both to the wing of an 
airplane or missile or to an all-movable, rectangular control surface 

with no gap. The curves of figure 23 start at values of Xc *? - * corre- 

c 

sponding to those for the wing alone at — = 0. As increases for 

pa pa 

constant pA there is a forward movement of the center of pressure 
because of the loss of lift due to interference which is mostly effective 
on the rear of the wing. For the lowest effective aspect ratio of 2 
there is about a 4-percent forward movement of the center of pressure 
due to interference in the extreme case. For large effective aspect 

ratios the forward movement is not nearly so large. As the value 

increases for constant pA there is an asymptotic approach of the 
of pressure back to the wing-alone value. 



center 


Drag 


The results of the lift calculations can be used directly to deter- 
mine the effect of interference on the wave drag at zero angle of attack 
of a symmetrical, double-wedge, rectangular wing mounted on a circular 
cylinder. In fact, it is not difficult to derive the result that this ■ 
drag coefficient based on the area of the exposed wing is given simply as 


PCp _ 3 
(t/c) 2 


' I?*) * ( 2fsA ’ f) ( 1 ' ^s)MPA»o*) 06) 


Equation ( 36 ) is for the maximum- thickness position at the midchord, 
although other maximum-thickness positions can easily be handled. 


The values of the drag parameter given by equation ( 36 ) have been 
calculated for effective aspect ratios of 2, 3* and 4 using the values 
of from figure 22. The results of the calculations are shown in 


figure 24 wherein 


(t/c) 2 


is plotted against 


c 

— — for constant values 
pa 


C bUn 1 

— = 0 — is equal to the wing-alone value of 4. 

Pa (t/c) 2 


of pA. At 
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As — increases, there is an increase in drag which in the extreme case 

pQ, 

amounts to about 4 percent for (3A = 2, and thereafter there is a steady- 
decrease in drag until at -p- = 12, the upper range of the calculations, 

p8. 

there is about a ^-percent decrease in drag from the wing-alone value. 

For other effective aspect ratios the changes are inversely proportional 
to the effective aspect ratio. It can thus be said for rectang ular 
wings of the present type with £A>2 and —-<12 that no important 

effects of interference on the mini mum drag occur. For missile-like con- 
figurations for which the wing contributes only a small part of the drag, 
the effects of interference on minimum drag will be negligible. For 
wing-body combinations having a large part of the wing lying b ehin d the 
Mach waves from the leading edge of the wing-body junctures, there can 
be appreciable effects of interference on ml nlrnimi wave drag. 


CONCLUDING REMARKS 


A theoretical method has been developed that makes it possible to 
obtain the pressure field acting on a wing-body combination composed of 
a wing with supersonic leading and trailing edges and a circular body. 

The method has been applied to the calculation of the pressure fields 
acting on combinations of bodies and rectangular wings with the bodies at 
zero angle of attack. The exact calculations are for combinations for 
which the effective aspect ratio of the exposed wing panels joined together 
is greater than 2, and for which the effective chord-radius ratio is 4 or 
less. For the family of combinations for which the effective chord-radius 
ratio is k, it was found that the blanketed area of the wing acts effec- 
tively at zero angle of attack and served primarily to support lift carried 
over from the wing onto the body. It was determined that for this family 
of combinations the body was somewhat less than 50 percent effective, com- 
pared to a perfect reflection plane, in reflecting pressure waves back 
onto the wing. 

By comparing the contributions to the span loading of the various 
Fourier components , it was determined that for rectangular wing-body 
combinations with effective chord-radius ratios of 4 or greater, the 
loss of lift due to interference is associated princip ally with the first 
component . On the basis of this fact it is possible to obtain simple 
asymptotic formulas for wing lift and wing center-of -pressure positions 
for large values of the effective chord-radius ratio. The results of 
the exact calculations and of those based on the asymptotic form ula s 
enabled the construction of design charts show ing the wing lif t and 
center-of -pressure position as a function of effective aspect ratio and 
effective chord-radius ratio. These charts are applicable to wings or 
control surfaces deflected in the presence of a round body. The char ts 
show that as a result of interference the lift of the wing panels in com- 
bination can be reduced as much as 15 percent from that of the half -wings 
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joined together, and that the center of pressure of the wing can move 
forward as much as 4 percent of the chord. Ho large effects of inter- 
ference on the wave drag of symmetrical double-wedge rectangular wings 
with the maximum thickness at the midchord were found. 


Ames Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Cal if ., January 8, 1952. 
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APPENDIX A 

DETERMINATION OF FUNDAMENTAL SOLUTIONS FOR FOURIER " 
COMPONENT COMBINATIONS 


Von Karman and Moore determine the potential for a body of revolu- 
tion by distributing along the axis of the body fundamental solutions of 
an axially symmetric type. This naturally suggests the use of fundamen- 
tal solutions varying a.3 cos 2n0 for the Fourier component combinations 
for which the potential varies as cos 2n0. Such fundamental solutions 
satisfy the boundary condition given by equation (5)» The bo undar y con- 
dition given by equation (6) can be satisfied by making -the fundamental 
solution imaginary in the forecone, that is, by m aking it proportional 

to [ (x-£ ) 3 -r 2 ] m / 3 where m is an odd integer. These conditions make 
the desired fundamental solutions of the form 

cos 2n0 [ (x-i i|r(x-£,r) 

where ^(x-i,r) is a function to be determined in such a manner that the 
fundamental solution obeys the wave equation. The parameter £ is the 
position on the x axis where the ■ fundamental solution is located. In 
the case of a wing-body combination, £ will vary from -1 to + 00 . How- 
ever, a value of 1=0 is used in the derivation of the fundamental 
solution so that it will be at the origin. If, as a simple possibility, 
i)r(x,r) is taken as a power of r (or x) and the fundamental solution is 
substituted into the wave equation in cylindrical coordinates (see fig. 1 
for coordinate system) 


*00 


^rr + r + r 2 *xx 0 


(Al) 


two possible sets of fundamental solutions are obtained 


cos 2n0 



r 2 n 


(A2) 


pn 

cos 2n0 r (A3) 

(x 2 -r 2 ) 2n+ 5 

The first solution is singular along the body axis (r = 0), and the 
second solution is singular along the Mach cone (x = r). Derivatives of 
the above solutions with respect to x are again possible fundamental 
solutions . 


A Bet of fundamental solutions having its singularities along the 
body axis would be better than a set having its singularities on the 
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Mach, cone since it was desired to get the external flow field about a 
combination. However, when the foregoing fundamental, solutions with the 
singularities on the axis were used, it was discovered that the functions 
defining the strengths of the distributions along the body axis had 
u n m an ageable singularities. In fact, only by using the (2n-l)'th deriv- 
ative with respect to x of these solutions was another set of funda- 
mental solutions obtained which gave manageable singularities in the 
strengths of the distributions along the axis. The set of fundamental 
solutions, <P 2n *(x,r,0), 


1WU,r,e) - (x 2 -* 2 ) 211 '* m 

are those adopted for this report. For n = 0, the partial differentia- 
tion convention of equation (Ak) is taken to indi cate an integration 
with regard to x. 
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APPENDIX B 

DETERMINATION OF AXIAL STRENGTH FUNCTIONS FOR FOURIER 
COMPONENT COMBINATIONS 


The distribution along the axis of the fundamental solution is deter- 
mined by the boundary condition of equation (9)* that is, the distribu- 
tion is to produce a given velocity amplitude function f 2n (x). If the 
distribution of cp 2n *(x,r,0) (equation (A4)) along the body axis 
from -1 to 00 to produce a step function of f 2n (x) were known, then 
the axial strength distribution to produce the given function f2n( x ) 
could be closely approximated by step functions as shown in figure 25 * 

Let h 2n (x-l) be the distribution along the axis from -1 to 00 to pro- 
duce a unit step in f 2n (x) at the origin. Then to obtain the axial 
strength function to produce the stepwise approximation to f an (x) shown 
in the figure, the following sum must be formed: 





If g 2n *(x-l) is the strength distribution along the axis for the actual 
f pp (x) distribution found by letting the Ax's approach zero, then 
since ^?n( 0 ) — 0 

®2n*( 

with the definitions 


x-i) - r ^2n^ x ”^-”^ )^Hn 1 ( 5 ) d£ 


(Bl) 


g2n*( x “ 1 ) = g2n( x ) 


(B2) 


and 


h 2n ( x -l) 


2 an (2ni ) 

(knl) 


^2n( x ) 


(B3) 


there is obtained 


g2n(z) =- — {Oil f" M^x-S) f ffl ' (£) d£ (B^) 

(4nl ) Jo 

The changes introduced by equations (B2) and (B3) are convenient for 
future work. By equation (B3) the function Mgn (x) is a simple fraction 
of hpn(x-l), the axi r 1 strength distribution necessary to produce a unit 
step in f a (x). The functions Myn (x) have been studied and evaluated 
numerically in- reference 5> and the numerical values are tabulated in 
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table I. A plot of these functions is shown in figure 5* These func- 
tions , together with equation (Bh) , give a simple means of obtaining the 
axial strength distribution by means of numerical or graphical integra- 
tion for any desired f 2n (x). The possibility of finding universal 
functions Mpn(x) for doing this is a consequence of the fact that the 
body is cylindrical. In the Karman -Moore method no such s tuple rela- 
tionship is possible because the body is of arbitrary shape. 
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APPENDIX C 

ALTERNATE METHOD OF OBTAINING BODY PRESSURE COEFFICIENT 
OF FOURIER COMPONENT COMBINATION 


A method of obta i n i ng the body pressure coefficient involving only 
one numerical integration is now given. Consider a disturbance at the 
origin formed by a bump on the surface of the body as shown in figure 26 . 

f (x) 

The nondimens ional amplitude of the bump — is to be zero except 

just at the origin where the area under the fpn( x ) curve is unity, that 
is, the disturbance is a ’’delta function". The amplitude is to vary 
as cos 2n0. Now if the pressure coefficient in the wing-body juncture 
Is -2 W 2 n (x) due to the unit disturbance at the origin, then the con- 
tribution to the pressure coefficient at a downstream distance x due 
to a distribution of disturbances from the origin to the point in ques- 
tion is 


5P„„ cos 

tx 


^2n(£ )^2n( x “£ ) ^ 


(Cl) 


This integral Is the effect of all disturbances forward of the point and 
Is additive to the effect of the disturbance exactly at the point. The 
local value is that given by the Ackeret theory, namely. 


P = 



(C2) 


where Is the local angle of attack 

f 2 n( x ) / % 

a = cos 2nd (C3) 

L V 

The pressure coefficient due to all disturbance is thus 

P 2n = | cos 2 n 0 f 2 n (x) - | cos 2 n 0 P f 2n (t JW^Cx-g) d| (CA-) 

V V J Q 

The advantage of this method of obtaining the pressure coefficient at the 
body surface is that the pressure coefficient is obtained directly as a 
function of the boundary condition, f^Cx), without first determining an 
axial strength function. 

The ¥ 2n (x) functions introduced to determine the body pressure dis- 
tributions are uni versal functions that have been studied in reference 5 . 
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They are presented in table II, and curves of them are presented in fig- 
ure 7. The values of W Q (x) have been tabulated by G. N. Ward in refer- 
ence 11. 

Several interesting facts are illustrated by figure 7* The ampli- 
tude of the pressure distribution due to a unit disturbance of the type 
shown in figure 26 has exactly 2n points of zero pressure, and these 
points are nearly evenly spaced. The large number of points of zero 
pressure for the higher harmonics makes their lift contributions much 
less significant than their pressure contributions. Another interesting 
feature of figure 7 is that the W^n(x) functions damp much more rapidly 
in the downstream direction for the higher harmonics than the lower har- 
monics. In fact, for n = 0 the damping to b percent of the maximum 
value occurs about 7 body radii downstream, whereas for n = 5 the same 
damping occurs in about 2 body radii. 
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TABLE I.- VALUES OF M~(x) 


x M 0 (x) 

M 2 (x) 

M 4 (x) 

M 6 (x) 

-CO 

-00 

-CO 

-CO 

.05 

: 

-531 

-.716 

.1 -1.321 

-.767 

.405 

1.334 

.15 

— 

— 

1.155 

.2 -.868 

-.166 

.902 

.735 

.25 

— 

— 

.279 

.3 

•131 

.831 

-.119 

.35 

— 

— 

-.417 

.40 -.538 

-304 

.560 

-.595 

.45 

— 

— 

-.675 

.50 

— 

.245 

-.655 

.6 -.378 

.450 

-.o4o 

-.417 

.7 

— 

-.253 

-.070 

.75 - — 

— 

— 

.095 

.8 -.292 

.457 

-.386 

.235 

.85 

— 

— 

.340 

.9 


-.434 

.4o6 

.95 - — 

— 

— 

-433 


.392 -.412 

-337 

.294 -.231 

-.115 

.188 .000 

' 095 

.088 .163 

.204 

.005 .211 

.197 

-.058 .164 

-.100 .069 

.020 

-.123 -.022 

.053 

-.126 -.073 

-.119 -.080 

-.071 

-.103 -.060 

-.038 

-.083 -.019 

-.002 

-.061 .011 

.021 

-.040 .027 

.027 

-.022 . 025 

.022 

-.007 . 015 


NOTE that Man fx - ) 




Wo(x) 

W 2 (x) 

w 4 (x) 

0.500 

0.500 

0.500 

— 


1 

.4319 

.7541 

1.644 


.3284 

.2887 -826 -706 -.940 -3 


.2550 

.706 

.2262 

-557 

.2014 

.402 

.1801 

-255 

.1616 

.128 

.1454 

.026 

.1313 

-.048 

.1189 

-.097 

.1080 

.0983 

.0898 

.0822 

.0754 

.0693 

.0638 

.0590 

-.122 

-.128 

-.122 

-.104 

-.082 

-.058 

-.035 

-.017 



















































Figure /. - Components of typical wing -body combination. 
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(a) fb) ' (c) 

Figure 2. - Decomposition of wing- body combination into component 

configurations. 



Figure 3.- Wing-body combination as the sum of wing atone plus 

interference combination. 





Figure 4.- Decomposition of interference combination into series of Fourier 

component interference combinations. 
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Figure 5. - Graphical presentation of M in (x) functions. 
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Figure 8. — Flow field produced by wing alone. 



Figure 9. - Variation of normal velocity induced at body surface 
by wing alone, unit body radius . 
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(a) x<i 



Figure U.- Singularities encountered in determining axial strength 

functions. 
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Figure 12 - Graphical presentation of axial strength functions. 
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Figure !4 - Singularities encountered in determining interference 

pressure coefficients. 
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(b) n=! 

Figure !5. - Interference pressure distributions of various Fourier components. 






x/fia -r/a+l 
(d) n*3 

Figure 15.- Concluded. 














Figure 17. -Pressure distribution on 


¥ 
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x/fia 

of combination of body and rectangular wing using four and six 
Fourier components. 
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Figure !9. - isometric drawing of pressure distribution acting on combination of body and 

rectangular wing . 
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Figure 20. - Span loading of various Fourier components 
acting on combination of body and rectangular wing 
having effective chord-radius ratio of 4. 
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Figure 2 1 . - Complete span loading for combination of body and 
rectangular wing having effective chord-radius ratio of 4. 





Figure 22. -Lift effectiveness for wing or control surface of combination of body and rectangular 

wing , 
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Figure 23.— Center of pressure for wing or control surface of combination of body and rectangular 
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Figure 24 - L ffect of interference on the minimum wave drag of a rectangular double - 
wedge wing mounted on a circular body, maximum thickness at wing midchord. 
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Figure 26. - Fourier component interference combination with 
"delta- function" protruberance at x=0. 
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